Shear Viscosity in a Gluon Gas 



OO 

o 
o 



(N 



> 

in 
d 

o 



Zhe Xu^ and Carsten Greiner^ 

'institut fiir Theoretische Physik, Johann Wolfgang Goethe- Universitat Frankfurt, 
Max-von-Laue-Str.l, D-60438 Frankfurt am Main, Germany 
(Dated: May 27, 2008) 

The relation of the shear viscosity coefficient to the recently introduced transport rate is derived 
within relativistic kinetic theory. We calculate the shear viscosity over entropy ratio r]/s for a 
gluon gas, which involves elastic gg — > gg perturbative QCD (PQCD) scatterings as well as inelastic 
99 999 PQCD bremsstrahlung. For Us = 0.3 we find rj/s = 0.13 and for Os = 0.6, ri/s = 0.076. 
The small rj/s values, which suggest strongly coupled systems, are due to the gluon bremsstrahlung 
incorporated. 



The elliptic flow measurements at the Relativistic 
Heavy Ion Collider (RHIC) indicate that the new mat- 
ter created is a nearly perfect quark gluon plasma fluid 
[H. Quarks and gluons should be strongly coupled [1]. 
The reason for it is still open. Attempts to understand 
the phenomena by using perturbative QCD (PQCD), 
which include elastic gg gg interactions, failed since 
gg — > gg interactions cannot even drive the system to- 
ward thermalization Q. Also, the shear viscosity to 
entropy ratio for elastic processes 0, 0| is much larger 
than the lower bound rj/s = I/Att found from strongly 
coupled supersymmetric Yang-Mills gauge theory using 
the AdS/CFT conjecture However, the importance 
of PQCD bremsstrahlung is raised in the "bottom-up" 
thermalization scenario 7] and the physics regarding jet- 
quenching ^ . Recent calculations within the Boltzmann 
approach of multiparton scatterings (BAMPS) 0, [lo| . 
which includes PQCD gg <-!■ ggg bremsstrahlung, demon- 
strated that the gg <-> ggg processes are a factor of 5 more 
efficient for thermal equilibration and pressure buildup 
than the gg — > gg scatterings, whereas their transition 
rates are almost the same. It raises again the question 
whether PQCD interactions can in fact explain that the 
quark gluon plasma behaves like a "strongly coupled" 
system with a small shear viscosity to entropy ratio. 

In this Letter we first derive a useful formula for the 
shear viscosity coefficient within relativistic kinetic the- 
ory and then calculate the shear viscosity to entropy ra- 
tio for a gluon gas in thermal equilibrium, which includes 
PQCD gg — > gg and gg <~+ ggg interactions. 

For a gas the phase space distribution of particles 
f{p,x) satisfies the Boltzmann equation 



(1) 



where — (1, p/E') and / denotes the collision term that 
determines the change in / due to interactions among 
particles. If the gas is locally in kinetic equilibrium but 
still away from chemical equilibrium, its phase space dis- 
tribution is described by the distribution 



where \{x) is the fugacity, I3{x) is the inverse of temper- 
ature T{x), and u^{x) denotes the four- velocity of the 
medium. The =F sign applies to bosons and fermions, 
respectively. 

The viscosity can be extracted from the response of 
the medium to a small deviation from /o 



/ = /o + /o(l±/o)/i 
by comparing the stress tensor defined as 
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with its expression in the Navier-Stokes approximation 
T"" = (e + P)u''u'' - Pg"'' + r^'' , (5) 

where 

1] denotes the shear viscosity and ^ the bulk viscosity. 
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u^u^, V" = A"^90 where g 



diag(l, —1, —1,-1). For ^ Landau's definition of u'^ is 
used llll . In the rest frame of the medium we obtain 
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In a conformal theory, which can be applied to the (mass- 
less) gluon gas at high temperature, T^^ -\- Tyy + T^z = 
Too — £ — 3P so that the bulk viscosity vanishes. 
With ^ and using 

d^fo = /o(l ± /o) [9,. In A - d^{(3u,pn] , (9) 
the left-hand side of the Boltzmann equation ([1]) becomes 
v'-dj = /o(l±/o)[l + (l±/o)/i±/o/i] 

+/o(l ± .fo)v^df,h 
^ fv^d^ In A - /o(l ± fo)v^d^ {(3u,pn , (10) 
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where the neglected terms are of second order in spa- 
tial gradients when assuming that fi in ([3]) only con- 
tains terms of first order in spatial gradients. Following 
di_iT'^'^ = 0, Eq. PH)) is rewritten in the rest frame 

v^'d^!^ >^a^lnA + /o(l±/o)|^ 

X [9t(lnA)i;e-|-3a^(lnA)P*r''*] 
+/o(l ± fo)PEv'v^ 

^{d,u^ + d,u') ~ y^jW ■ u , (11) 

where i,j = 1,2,3 and the equation of state e = 3P 
is used for massless particles. For A = 1 Eq. pTI) is 
identical with the term derived by Arnold et al. in [4| . 
Integrating Eq. ([1]) over momentum using (jlip gives 
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(2^)3^"^/^/ = J"5,(lnA), (12) 



where 
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is the particle density. In the derivation for ([T^ we as- 
sumed that there is no particle flow in the rest frame 
following Eckard's definition of ll|. Integrating ([1]) 
by weight gives 
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where (•) denotes the average over particles. We then 
obtain the relation 



where 



J2R'" + jndtilnX) 



(15) 



= ^''^ 

is the total transport collision rate, which was introduced 
in [lol | as the characteristic quantity describing momen- 
tum isotropization. In kinetic equilibrium, Eq. (|15p be- 
comes an exact equation. 

Inserting ([TS]) into ([7]) we obtain 
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J2R''^ + indt{lnX) 
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This expression constitutes our major formula calculating 
the shear viscosity coefficient ry and gives a direct corre- 
spondence of T] to the transport rate i?'' : rj is inversely 



proportional to the sum of the total transport collision 
rate and the chemical equilibration rate, and is roughly 
proportional to the energy density. If the chemical equi- 
libration is governed by 2 3 processes, 



-ndt{lnX} = 3 



d'p 3 
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which might even become negative for systems with over- 
saturation. However, if the system is not far away from 
chemical equilibrium, the total transport collision rate is 
most dominant and, thus, determines the shear viscosity. 

If / or /i in ([3]) is known and /i ^ 0, ry can be in 
principle calculated using ([7]) . A way to get /i is to solve 
the linearized Boltzmann equation 

fv^'d^ In A - /o(l ± fo)v^'d^{l3u,p'') = /(/i) (19) 

as a variational problem for /i [llj . Often used is a sim- 
ple ansatz 0, [3 fo'^ A such as the function (PTjl below. 
One notices that 77 obtained from ([7]) using / from is 
identical with that from (|17p . Another more complicated 
method to get / is to solve the Boltzmann equation ([1]) 
numerically performing extensive transport simulations, 
which is, in principle, more reliable than the first method 
except for numerical uncertainties. Since the spatial gra- 
dients are needed in ([7]) and their extractions from trans- 
port simulations are difficult, it is at present more con- 
venient to use (fT7|) to calculate r]. Such calculations will 
be presented in a forthcoming publication. 

In this Letter we consider a static particle system, 
which is initially disturbed from equilibrium and is relax- 
ing again to thermal equilibrium. For a Boltzmann gas 
in chemical equilibrium / is assumed to have the form 

M 



in the rest frame. For small x 



f 



(20) 



(21) 



Comparing (PT|) with ^ fi = ~ f /3 ^ , which satisfies 
Eq. with appropriate x- Using (pTj) in ([TT]) we cal- 
culate r] in thermal equilibrium as the limit when letting 
X — > 0. In this limit Eq. (jl7[) becomes an exact equation 
identical to Eq. ([7]). We obtain 
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where Z*'' is the inverse of the total transport collision 
rate and is called the mean transport path. 

We now apply Eq. (|22|) to calculate the shear viscosity 
in a gluon gas. For the sake of simplicity we assume glu- 
ons to be Boltzmann particles. The boson enhancement 
is neglected, which is a good approximation for gluons 
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at high temperatures. Interactions among gluons include 
elastic gg gg PQCD scatterings and inelastic gg <-^- ggg 
PQCD bremsstrahlung processes, which are screened by 
the Debye mass 

/(fip 11 1 
TTr-y^-^cfo ^ TT '^G Nctts— , (23) 

where da — 16 denotes the gluon degeneracy factor and 
Nc — 3. The Landau-Pomeranchuk-Migdal suppression 
of bremsstrahlung is taken into account as a lower cutoff 
in the momentum of the radiated gluon [3| • The matrix 
elements of the transitions can be found in [l3| • We 
obtain ^ i?'' in thermal equilibrium as the term in 0th 
order of x Eq. (fT6|) when using Eq. (|2T|) for the 
integrals. Because the collision term is additive, J2 — 

Dtr I Dtr _|_ Dtr whprp — /?^^ 

^gg^gg^^gg^ggg^^gggy^ ^ ^ggg^gg ~ 3 gg~*ggg 
due to detailed balance [1(|. 

The shear viscosity to entropy ratio at equilibrium is 

^ = (5/3 E^") ' = {^^^I'^-^g + f K™) ' . 

(24) 

where s = |/3e is used. Because the transport collision 
rates scale with the temperature, (iR^gg^gg, PR%->ggg 
and, thus, rj/s depend only on the coupling constant a^- 
The upper panel of Fig. [T] shows the shear viscosity to 
entropy ratio for gg — > gg and gg -s-^- ggg processes, re- 
spectively, rj/s for gg <-!■ ggg processes is roughly a factor 
of 7 smaller than that for the elastic collisions, which im- 
plies that compared to the elastic scatterings the PQCD 



inspired bremsstrahlung is the leading process in relax- 
ing the system to equilibrium. For as = 0.3, which 
might be appropriate at RHIC energy, rj/s = 1.03 for 
gg —^ gg only and rj/s = 0.13 when including gg ggg 
processes. To match the lower bound oi -q/s = l/47r 
from the AdS/CFT conjecture 0] c«s = 0.6 has to be 
chosen. The -q/s ratios in the two cases correspond to 
the mean transport path Z''' = 1/^i?*'' = 0.32 fm and 
l^' = 0.2 fm at T = 1//3 = 400 MeV for as = 0.3 
and as = 0.6, respectively. From the collision rates 
shown in the lower panel of Fig. [T]with n = doT^ /t^^ 
we obtain (vi-c\(t gg^gg) = Rgg^gg/n = 0.82 mb and 
{vi-c\o-gg^ggg) = Rgg^ggg/u = 0.57 mb for as = 0.3, and 
(I'loiCTgg^gg) = 1.27 mb aud (vreiCTgg^ggg) = 0.73 mb for 
as = 0.6 at r = 400 MeV. Hence, perturbative interac- 
tions can drive gluons to behave like a strongly coupled 
system with a small rj/s ratio at RHIC. 

From Fig. [T] we also obtain that Rfg^gg/Rgg^gg = 
0.36(0.46) and Rl'g^ggg/ Rgg-,ggg = 2.1(2.7) for a, = 
0.3(0.6), respectively. The wide difference in the behav- 
ior of the R^'' / R ratio for the gg — > gg and gg — > ggg 
processes is essential for the different contributions to 
77/s. Because i?'' contains an indirect relationship with 
the distribution of the collision angle 0, we decompose 
the transport collision rate to 

i?-' = Ai 7i(wrciCTf ) , i ^ gg gg , gg ^ ggg , (25) 

with = J dcji siv? 9 defined as the transport cross 
section [3| and Ai being a multiplication factor. Fig- 
ure [2] shows Pn{v^ciO'l'^) and Ai as function of a^. For 
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FIG. 1: Upper panel: Shear viscosity to entropy ratio for 
gg gg and gg ^ ggg processes, rj/s for gg gg is divided 
by a factor of 7. Lower panel: Collision rate to temperature 
ratio for gg gg and gg ggg processes. 
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FIG. 2: Scaled transport cross section (upper panel) and 
transport collision rate to scaled transport cross section ratio 
(lower panel) for gg —^ gg and gg — > ggg processes. 
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,)/(«,eia<,3^g<,) = 0.3(0.35) 
0.71(0.78), which in- 
values gg ~^ gg favors 
-^ ggg favors large-angle 



and {vrcicrl^g^ggg )/{v rclO'gg^ggg) 
dicate that for the chosen 
small-angle scattering and gg 
bremsstrahlung (loj . 

The factors Ai in Eq. ([25]) have weak dependences 
on as and are around the isotropic distribution val- 
ues. (For isotropic angular distribution Agg-^gg =9/8 

and Ann-^nnn 



tain Aa 



g = 27/16.) For as = 0.3(0.6) we ob- 
j/Agg^gg = 2.5(2.7), whlch are significantly 



is fitted by (0.68 



larger than 1. 

Parametrically, PR^gg^gg 
2.8as) Qf^ (In Qfg)^ from as = 0.001 up to 0.3, whereas 
(3Rgg^gg ~ (3 — 36as/TT)as for as < 0.1. Thus, 
Rfg_gg/Rgg-tgg ~ ©(ctfi) for Q!s < 0.1, whlch indicates 
again small-angle gg — > gg scatterings, because for those 
scatterings R^^g^gg/ Rgg^gg ~ 0'^ with 0'^ « W I ~ 
'im'jy/sm ''^ as where q is the momentum transfer. 

Including the Landau-Pomeranchuk-Migdal effect for 
the PQCD inspired bremsstrahlung within the Bethe- 
Heitler regime (iRgg^ggg ^ ^^(Inas)^ 9]. For small- 
angle bremsstrahlung PR%^ggg should be of the order 
a'lihias)^ , which is clearly not the case for as > 0.01 as 
seen from Fig. [2] as within this as range the transport 
cross sections for gg — > gg and gg — *■ ggg have almost 
the same as dependence, namely ©[a^ (Incts)^]. Thus, 
Rfg^ggg/Rgg^ggg ^ for as > 0.01, which shows 

that in the chosen as interval the bremsstrahlung fa- 
vors isotropic angular distribution. For smaller as, how- 
ever, the collision angle tends to be distributed forward. 
This can be observed in Fig. [2] where f]n{v^ciO'fg 
becomes steeper than I3n{vre\(y^gg^gg) at as < 
This implies that the "collinear" bremsstrahlung where 
R^gg^ggg / Rgg^ggg ~ 0{as) will occur at extreme small 
as and only then will have smaller contribution to the 
transport coefficients than the gg gg processes 

Finally, the ratio of the coUisional width to 



0.01. 



16|. 

the 



mean 

Raaa— 



gluon energy, r/{E) = {Rg 



Ro 



,)/3/3, can be calculated from the lower panel of 



1.5i?„ 



due to detailed bal- 



Fig. [T] where Rggg^ 
ance. For as = 0.3(0.6) we obtain r/{E) = 0.5(0.69). 
These ratios are smaller than, but, close to 1, which indi- 
cates that PQCD with as = 0.3 — 0.6 is at the edge of its 
applicability. For larger values of as one obtains F > (E) 
and the PQCD calculations with on shell kinematics are 
no longer applicable, as for such a strong coupling regime 
the Heisenberg uncertainty principle has to be taken care 
of by a full quantum transport treatment 



in 



The higher order processes such as ggg ggg and 
gg <-> gggg will certainly modify the total transition rate. 
However, their contributions are suppressed by higher or- 
der of as [isj . Because the full diagrammatic many-body 
theory for higher order collisions becomes rather complex 
[3], at present, the incorporation of higher order multi- 
particle interactions is considered in a phenomenological 
manner |20|. 



In summary, the shear viscosity is derived within rela- 
tivistic kinetic theory, which is proportional to energy 
density and inversely proportional to the total trans- 
port collision rate. We calculated the shear viscosity 
to entropy ratio rj/s for a gluon gas, and found 77/ s — 
0.13(0.076) for a, = 0.3(0.6). Perturbative QCD inter- 
actions can drive the gluon matter to a strongly coupled 
system with an 77/s ratio as small as the lower bound from 
the AdS/CFT conjecture. The PQCD inspired gluon 
bremsstrahlung is responsible for small 77/s ratios and, 
thus, can explain that the quark gluon plasma created at 
RHIC behaves like a nearly perfect fluid. 

Z. X. would like to thank N. Su for fruitful discussions. 
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